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ABSTRACT 


It is suggested that the distribution function approach be used for 
describing the turbulent chemical reaction in order to avoid violation of 
an invariance property noted by O’Brien: for zero molecular diffusion in 
homogeneous turbulence, single point statistical functions of the concen- 
tration field are independent of the motion. This property and a similar 
’’white noise” invariance for two point statistical functions are shown to 
apply to arbitrary isochoric motions when the initial concentration field 
is statistically homogeneous with respect to position. Finally, it is ob- 
served that the use of a reciprocal concentration avoids the problem of 
preserving the invariance in turbulence closure schemes. 
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I. INTRODUCTION 

0*Brien 1,2 has recently noted an invariance property of the model second 
order irreversible chemical reaction in turbulent flow 3 

ft (x,t>r(x, t) =t-cr 2 (x, t) (i) 

where ft (x, t ) 5 { B/Bt + U(x, t)«V - DV 2 }, F is the mass concentration in 
either mass or molar units (with appropriate dimensions of C), U (x, t ) is the 
incompressible mass-average velocity field, D is the constant molecular dif- 
fusivity of the species F with respect to the mixture, and C is a constant reac- 
tion rate constant. The restrictions that U be unaffected by F or by temperature 
changes and that D and C be constant usually require low F and either small 
heats of reaction or large heat capacities. That no reverse reaction or creation 
term appears in (1) usually restricts one to the initial stages of reaction, unless 
T is extremely dilute in the product of reaction. We adopt the outlook that F 
is a passive contaminant. 

The invariance property in homogeneous turbulence is that for D = 0 all 
statistical functions <(r n (x , t )^> n = 1, 2, . . . are independent of the function U. 
Here the angular brackets denote ensemble average over all realizations of U . 

In homogeneous turbulence we require (in this paper) that statistical functions of 
U and F are independent of position x . In the more general case, when U is 
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an arbitrary isochoric motion with no flow across the boundaries, the invariance 
requires that for D * 0 all volume averages of P^x, t ) be independent of U. 

We will not discuss the latter case here. 

One can show the invariance directly from (1) by setting up an unclosed 
hierarchy of statistical moment equations in homogeneous turbulence (or of 
volume averages in the general case). This hierarchy is independent of U . 
Alternatively, in Lagrangian coordinates the fluid motion enters (1) only through 
the Lagrangian form of the diffusion term. For D = 0 the behavior of P on any 
trajectory or particle path is independent of U , and the invariance follows from 
the requirement that P is initially statistically homogeneous, the identity of 
Lagrangian and Eulerian volume averages in isochoric motion, and an ergodic 
assumption on the identity of statistical and volume averages. Later, this last 
assumption will be relaxed. 

The purpose of this paper is to see how the invariance may be preserved 
in closure theories. We will present arguments for the distribution function 
approach and look at the one and two point probability density functions of P for 
D = 0. Then we will show that a simple change of variables can remove the 
problem. 

II. PRESERVING THE INVARIANCE 

There are two types of dynamic nonlinearity of (1) leading to closure diffi- 
culty of the moment equations. First, the nonlinearity of the reaction term and 
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statistical initial and boundary conditions produce coupling between moments of 
all orders of r . This difficulty is formally the same as that of the turbulence 
problem itself (closure of the Navier-Stokes equations), except that the coupling 
in (1) is local in x -space. Second, randomness is induced by the turbulence, 
analogous to the nonreactive passive scalar problem. Of course, the nonlineari- 
ties are not independent, as the coupling in the reaction term will be induced by 
turbulence alone if the initial and boundary conditions of r are certain. 

Many of the analytical methods in turbulence and turbulent convection theory 
have been treated as procedures for forming mom <nt closures regardless of the 
nature of the dynamic nonlinearity. The direct interaction (DI) approximation 4 
and the Lagrangian history direct interaction (LHDI) approximation, 5 as they 
have been applied to (l), 2 * 6,7 clearly fall in this category, as do the more clas- 
sical approaches of iteration expansions and cumulant discards. 1,8 The phe- 
nomenological theories have not been considered for this problem, since the re- 
action term is in no sense a transport term, but rather is purely dissipative. 

If the moment closure does not explicitly involve the velocity field in the 
reaction term, such as an ad hoc cumulant discard hypothesis applied to moments 
of T, the zero diffusion invariance will not be violated, although the decay rates 
of the moments of r and the influence of the velocity field for D ^ 0 may be 
given incorrectly. The DI approximation applied to (l) 6 does involve the 
velocity field in the reaction term but fails to preserve the invariance. 2 The 
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reaction is not treated exactly, as the model for which the DI approximation is 
the exact solution in the reaction problem involves artificial dynamical coupling 
In both the convection and reaction terms. However, the LHDI approximation 
(and the abridged version) does preserve the invariance. 7 Its success can be 
traced to the fact that the averaged response function G for infinites simal con- 
centration perturbations retains two important properties of the Lagrangian form 
of (1) for D = 0: first, all fluid particle trajectories are independent ( G is non- 
zero only for perturbation and response on the same trajectory), and second, all 
trajectories behave identically (G is independent of labelling time) so that G is 
independent of the velocity field. Despite this success, there is no guarantee of 
the accuracy of the LHDI's treatment of the reaction term. Furthermore, because 
of the presence of Lagrangian velocity correlations, we regard the LHDI to be 
too difficult in practice. 

If one takes the point of view 9 that the turbulence, acting through the random 
coefficient U, should provide the only stochastic indeterminacy in (1), and that 
the effect of the reaction should be treated exactly, as the diffusion term is 
customarily done, then any such closure method should retain the zero diffusivity 
invariance. Actually, the reaction need be treated exactly only in the limit D = 0. 

It is clear that if the reaction term is to be treated exactly, the moment 
closure approach is not appropriate to the study of (1). For example, any 
iteration procedure based upon treating the turbulent velocity as a perturbation 
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is very difficult in that it requires solution of a nonlinear and nonhomogeneous 
partial differential equation at each level of iteration. Even for D » 0 one is 
confronted with a sequence of Riccati equations for P . 

On the other hand, the distribution function approach seems appropriate 
because the reaction term causes no extra conceptual difficulty above that of 
the nonreactive passive scalar problem. The closure problem here is not due 
to the nonlinearity of the dynamical equation, but to nonlocalness of the dynamical 
variables. The reaction term is completely local in x -space, whereas the dif- 
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fusion and convection terms are not. 

In the next sections we briefly outline the distribution function approach and 
examine the zero diffusion invariance in greater detail. 

III. THE DISTRIBUTION FUNCTION APPROACH 

Wo now sketch out the distribution function approach to justify our preceding 
remarks and to focus on its closure difficulty. The overall objective is to obtain 
at time t the one and two point probability density functions (pdf*s) of P defined 
by 


p* (0) d ® p {0 < r (x, t) < 0 + d \py 


( 2 ) 


P *1 ,* 2 dl Ai #2 = 1 ® % - r (*r t) , <p a < r(* 2 , t) < i /» 2 + dip 2 > 
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where P denotes probability. These pdPs will allow us to calculate mean decay 
rates and two point correlations of r and may be>obtained by reduction from the 
Ml distribution functions. Let P* [0 (x) | u (y) ] , or P fc [if/ |u] , denote the con- 
ditional probability functional at time t for the concentration field F, given the 
velocity field U, i.e. 

P* [«//(x) | u (y) ] 8^(x) * P (0(x) $ F(x, t) < \jj (x) + d^(x), V x i 

u (y) S U(y, t) < u(y) + du(y), V y) 

where the conditional events are to the right of the vertical bar, and 8 if/ is ele- 
mental volume in f-space. Strictly speaking, the probability density functional 
does not exist in infinite dimensional phase space where 8 if/ cannot be defined, 
and must be interpreted as the joint pdf of T at an arbitrarily large number of 
points x k . A Liouville or continuity equation for P* [if/ |u] may be derived from 
(1) with the method used by Edwards 10 ; we get 



[Ct/' 2 (x) + u (x) * Vi//(x) -DV 2 t//(x)] P^if/M (3) 


where 8 / 8 denotes functional differentiation. 

To obtain an equation for P* [if/] , the marginal functional distribution, multiply 
(3) by the full pdf of U for all time, i.e. by F M , defined by 
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F [v] $v » P {v (y, s) i U (y, s) < v (y, s) + dv(y, s), Vy and 0 $ Vs < T) 

where T is some finite time greater than any considered in the problem, and 
then perform a functional integration over all v * The result is similar to (3) 
except for an integral over the convection term involving the joint distribution 
of ij/ with u. The ability to evaluate this term at any t is tantamount to 
knowing F O] , whereas in practice we can only deal with a few moments of , 
The procedure of Ref. 9 is to expand P* [y |u] in the convection term about 
P* [v|/] before carrying out the v integration, in order to obtain an expression 
for ? t [y] in terms of moments of v. 

» 

The single point pdf is obtained by multiplying (3) by F M and integrating 
over all v- and over \\ / for all points except x, The result is 



P* ( \p ) ” - D lim 
x' x 



dxp' 0' P* 


,x' 


OA 00 


/ d ' u 


(4) 


where the convection term, which involves the joint distribution of y with u , 
must be approximated, say as in the preceding paragraph. For the particular 
case of joint statistical homogeneity of f and U, i.e. P*,* (0> u) is independent 
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of x» v the convection term in (4) is zero. If D were also zero, then P*^) may be 
calculated from the result for zero motion in accordance with O'Brien's invariance. 
We shall discuss this in more detail later. Unfortunately , there is also a closure 
problem associated with the nonlocalness of the diffusion term. An unclosed 
hierarchy for these multipoint pdf's can be set up, with each equation similar 
in form to (4), and with the reaction term treated exactly, as m the left hand side 
of (4). For example the two point pdf obeys the equation 





u u * 


WvVv u ) + v x 


W'r'/V “)] 


Higher order equations may be obtained by induction. This sequence is analogous 
to that for the ordinary turbulence problem. 11,12 

In the ordinary passive scalar problem, for which C = 0, the nonlocalness of 
the diffusion term could be removed by using the spatial Fourier transform of 
(1). BJc for C / 0 the reaction term in wavenumber coordinates is a convolution, 
and closure on that term may lead to violation of the D= 0 invariance. In 
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^-coordinates a closure must be made on the diffusion term, but the reaction 
may be treated exactly and will preserve the invariance. However, for small 
C, large D, or in regions of phase space where |V<i/| is large, the invariance may 
not be as important, and in those instances it is more important that the diffusion 
be treated exactly. 

We suggest the following method for dealing with the diffusion term. First, 
for the velocity term, make a suitable closure forP* [y \ u] . Then, base the 
reduction from P l [y] for use in the diffusion term, say at the two point level, 
upon the formulation of (1) in k -space, where the diffusion may be treated exactly 
at the expense of exact treatment of the reaction. This approach would guarantee 
correct handling of the reaction term in (4) when diffusion is weak, with the dif- 
fusion term being treated properly when diffusion is dominant. 

IV. THE SINGLE POINT PROBABILITY DENSITY 
FUNCTION FOR D = 0 

Let us now consider the single point probability density function (pdf) of 
r in the limit D = 0 to see what form must be retained in any distribution 
function approach that satisfies O'Brien's invariance. In this section, the term 
homogeneity means that the attribute in question, say r (x , t ) , is statistically 
homogeneous with respect to position: the probability density function of the 
attribute is independent of the position of the attribute. 
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First, examine the case when there is no motion, U(x, t) = 0. Then (1) reduces 
to 

^■r(x, t> +cr ! (*,t) = o (5) 

which has already been studied by O’Brien. 1 The solution of (5) is 

r(x,t>»2‘ r(x,o> = r(x,o> [x 4 ct r(x,o>]-‘. 

> , . 

The single point pdf of F (x , t ) can be related to the initial distribution of F by 
the identity 


( 6 ) 

Now introduce the condition that T is initially homogeneous, i.e. P° is independent 
of x . From (6) the same condition prevails for F at time t , whereby 


n w = p x° (2-* *> 


P ‘ (I/-) = (1 - Ct ^)- 2 P° (—l—). (7) 

We now observe that P fc is the same as (7) for D = 0 in homogeneous tur- 
bulence if P° is independent of x. This is clear from (4). Alternatively, the 
single point moment equations, which are independent of U in accordance with 


10 


12 


O f Brien v s invariance, are the same as those for (5); and because the pdf, or 
equivalently the characteristic function, is determined completely by the moments, 
the equality of the P's follows* 

Let us examine P* from another point of view* Rewrite (1) in Lagrangian 
form as 


- 3 - r« , t) + c r 2 (5, t) = D c-‘ (5, t) : V 5 r (5, t) (8) 

where § is position in the Lagrangian or material frame, C _1 is the inverse Green 
deformation tensor 


C- 1 <C,t) = [(Vf- 1 ) • (7f- l ) tr lx,f(§,,)- 

the motion is given by x = f (§ , t ) with inverse § = f _l (x, t) and velocity 
U (x , t ) = Bf (f , t )/dt , the Eulerian and Lagrangian concentrations are related 
as 


r(x, t) = r<f«,t) f t)»r<e,t), 

and V and Vg are the Eulerian and Lagrangian del operators. 13 In analogy with 
(2) we introduce the Lagrangian single particle pdf of T 

P| 0//) d\p ~ P {\fj < r (§, t) < \p + di//}, 
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also the conditional Eulerian pdf of T for a particular fluid particle 

P x|g WO P {'A - P(x, t) < «// + di/f | g < f " 1 (x, t) < g +dg} 
and the pdf for the initialization of trajectories that pass through x at the time t 

F* (g) d 3 g s P {g < f “ l (x, t) < g + dg} . 

For the case D = 0 (8) has the same form as (5), i.e. 

ir( 5 ,t) + cr 2 ({,t) = o, O) 

and for homogeneous initial conditions P| = P ® = P° independent of x or § t 
where we have noted f (f , 0) = g . From (6) we obtain the result analogous to (7) 

p| (</.) = (l - ct 4>r 2 P° ( Yz^tj) = 0/0 (10) 

independent of § . The Eulerian pdf is recovered by summing the conditional pdf 
over all particles, 


P* 

X 




pt 

X 


('/') F* (g) d V 


( 11 ) 
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Note that P*| f ( 0 ) and ( 0 ) are pdf's for the same fluid particle g at the 
same time; they are therefore equal and, from (10) # independent of g. Conse 
quently, g integrates out of (11) with the result 


K <#) s9t M 

independent of x and of the motion, or 

P* (!P) = M = (1 - Ct ^)- 2 P° / Y^j j • (12) 


Recall the weakness of the assumptions used in obtaining (12) , Of course 
D = 0, C = constant, and P° ( 0 ) must be independent of x . The motion must also 
be isochoric so that (9) does not contain a dilatation term. However, we have 
made no assumption about the turbulence itself, such as spatial homogeneity 
of the velocity fluid. The derivation does not even require the motion to be 
turbulent, as (12) holds for each realization of the motion. We have had to 
introduce, however, an ensemble of initial P fields in addition to the ensemble 
of velocity fields. 


The following physical picture is presented by (12). Because P° is the same 


» 


for all particles, any particle that passes x at time t had the same P° in each 
realization of U and also the same age t. Because P * depends only on age for 
D = 0, it will be the same at x in each realization regardless of the motion, and 


the same for all points x . 
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Our conclusion for D = 0 is, then, that if at some time r becomes homogeneous, 
it will remain so independent of the subsequent motion, and P t may be calculated 
from the result (7) for zero motion. It would also be desirable to determine 
the circumstances under which P becomes homogeneous, although we will not 
investigate this here. (We would like to use the continuity equation for P t [i|/|u] 
to show that P tends toward homogeneity if the U field is completely homogeneous.) 

The Lagrangian formulation had to be used here because there is no counter- 
part to (6) in the Eulerian formulation. The order of the reaction is unimportant; 
in other instances the equivalent of (9) would be that P(§, t) be independent 
of other fluid particles so that (6) can be used. For no reaction we obtain the 
familiar passive scalar result P t (i //) = P° («//), i.e. all moments of P are constants 
of the motion. 


V. A "WHITE NOISE" INVARIANCE FOR D = 0 
In addition to the property of the single point pdf given by (12) for D = 0, 
there is a similar invariance for the two point pdf of F : if P is initially homogeneous 
and pairwise independent for all points, then it will remain so regardless of the 
subsequent motion. 


The proof is similar to that for (12), except for the use of P* x (i p v \ p 2 ), 
the two point pdf of P; Pi - (0, , 0, )» the two particle pdf of T ; P v fc \ 

x l* x 2'"l '*2 

(0i , 0 2 ), the conditional Euleriian pdf of r for two particles; and F^ tK (g x > g 2 )» 
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the joint pdf of particles passing through x l and x 2 at time t - all defined by 
analogy to the single point functions. The formula equivalent to (11) is 


P* 


,x. 




JT* 




i'Pv't’i) F * 1 ,x 2 (fti-8 2 )d 3 B, d 3 B 2 


JJ 




('/'[I 02 ) ^,,*2 (*1-82) d3 Bl d3 B 2 


(13) 


and that equivalent to (6) is 




(3 _t 0,. 3-* 0 2 ) 


3(3‘ l 0,, 2 -t i/< 2 
3 0/' l .</' 2 ) 


(14) 


where the Jacobian in (14) reduces to 


1(1 "Ct </-,) (1 -Ct ^)]- 2 


Under the conditions that P is initially homogeneous and pairwise independent 
for all points, 


,§ 2 (^^ 2 ) =P° (0 X ) P° (0 2 ) 

and (14) becomes 
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«V0 2 ) = ?' (0.) (0 2 ) (15) 

where P t is given by (10), With the substitution of (15) into (13), F* „ inte- 

* I, i* 2 

grates out and we obtain, using (12) 

p * t ,,* J (0i-0 2 )= :pt (0.) pt (02) (16) 

independent of the motion and of the points x t and % 2 . Consequently, homogeneity 
and pairwise independence is preserved. Results analogous to (16) also hold for 
multipoint pdf’s of any order. 

This invariance may be called a ’’white noise" invariance, as it implies for 
D = 0 that if the truncated power spectrum of P (in k -space) is initially flat, 
it will remain so, with all spectral components decaying at the same rate. Hence, 
once spatial correlations are lost, in the absence of diffusion they may not be 
regained. It would also be desirable, although we will not attempt this here, to 
be able to demonstrate with the continuity equation for P t [y |u] that P tends 
towards this equipartition. 

The DI approximation evidently does not satisfy this invariance, at least 
one cannot tell by examination of Eqs. (1,11) and (1,12) of Hef. 6. However, 
the LHDI approximation (and also the abridged version) does preserve it, as 
one can tell from Ref. 7 by using Eqs. (38), (39), (45) - (47), the properties of 
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G when D ** 0, and by noting that a property similar to Eq. (62) holds for the 
Fourier transform of the covariance of the generalized concentration field 
(normalized by (2tt)~ z x volume of wavenumber space). 

VI. THE USE OF SPECIFIC VOLUME 
Let us turn our attention to the reciprocal of the concentration of the con- 
taminant P, its specific volume 

t(x, t)*r - 1 (x,t). (17) 

Any moment of T depends on all the moments of r, and vice versa. From (1) 
the behavior of T is given by 


n(x, t) Y(x,t> =C~Y~^r |VY(x,t)| 2 . (18) 

The quadratic nonlinearity of the reaction term in (1) has been replaced at 
this level in (18) by a logarithmic nonlinearity in the dissipation term. Such a 
nonlinearity is actually present in the single point hierarchy equations for 
either (1) or (18), as we have 

Hr n = -ncr« + l - (n - 1) D (V log T) • (V T") 

HT n = nCY"- 1 - (n + 1) D (V log T) • (V Y"). 
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The standard practice (nonreactive passive scalar problem) for avoiding this 
dissipative term is to consider only the P hierarchy for multiple point correla- 
tions. After closure, certain points are made to coincide before substitution 
into lower order equations. The same approach could be used for C ^ 0, but 
one must be careful with the role played by the velocity field in the multipoint 
moments so that the O'Brien invariance is retained when the points coincide and 
D 0. Unfortunately, the log term appears in the first equation of the Y hierarchy. 
It is not apparent how to avoid dealing with it. 


Because the dissipation term in (18) disappears when D = 0, whatever type 
of statistical approximation we apply to this term cannot destroy O'Brien's 
invariance. For D = 0 in homogeneous turbulence, (1) and (18) give the following 
hierarchies of moment equations 

®l = - C © 2 <t>,=C 

» i i • • • 

and 

0„ :-nC0 ni1 <b ~ n C . 

n n + 1 n n - i 

9 9 9 9 9 9 

where @ n * <T n )> and \ 53 <(Y n ^> are averages defined over ensembles of initial 
data and velocity fields. Each hierarchy preserves O'Brien's invariance. How- 
ever, the $ hierarchy has the advantage of being closed at each level, and trunca- 
tion does not alter from the true values the growth of specific volume or its 
moments, contrary to the 0 hierarchy. Another advantage of dealing with (18) 
is that the conditions under which the model (1) is at all applicable, dilute liquid 
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phase reactions for which D is extremely small and Y is large, may allow 
the dissipation term to be linearized. 

A substitution similar to (17) can be used for irreversible reactions of 
arbitrary order N 4 1, Y * F l * N , yielding an equation equivalent to (18) for the 
(N - l)st power of the specific volume, 

vn. CONCLUSIONS 

We have suggested two approaches that could be followed to guarantee pres- 
ervation of O'BrienVs invariance and the correct rates of decay of single point 
moments when D « 0. Description by the distribution function approach is difficult, 
and we have not pursued it far enough to demonstrate practical use. The use of 
specific volume in place of concentration in the moment approach is simpler, but 
the experimentalist would prefer to see results in terms of concentration, and the 
singularity in Y for large times may be difficult to handle. 

The generalizations (12) and (16) of the zero diffusivity invariance for the 
one and two point pdf*s of V should, in retrospect, seem obvious on physical 
grounds. These properties should provide tests for the adequacies of either 
moment or distribution function approaches to either (1) or (18). 14 
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